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TECHNICAL MEMORANDUM NO. 1111 


THE. MONOPLANE AS A LIFTING VORTEX SURFACE 1 
By Hermann Blenk 


In Prandtl’s airfoil theory the monoplane was replaced "by a single 
lifting vortex line (reference l) and yielded fairly practical results. 
However, the theory remained restricted to the straight wing. Yawed 
wings and those curved in flight direction could not he computed with 
this first approximation; for these the chordwise lift distribution must 
he taken into consideration. For the two- dime nsional problem the transi- 
tion from the lifting line to the lifting surface has been explained by 
Birnbaum (reference 2) . In the present report the transition to the 
three-dimensional problem is undertaken. 

The first fundamental problem involves the prediction of flow, .pro- 
file, and drag for prescribed circulation distribution on the sti'aight 
rectangular wing (fig. 1(a)), the yawed wing fox- lateral boundaries par- 
allel to the direction of flight (fig. 1(b)), the swept-back wing (fig. 
1(c)), and the rectangular wing in slipping (fig. 1(d)), with the neces- 
sary series developments for carrying through the calculations, the 
practical I'ange of convergence of which does not comprise the wing tips 
or the break point of the swept-back wing. The second problem concerns 
the calculation of the circulation distribution with given profile for a 
slipping rectangular monoplane with flat profile and aspect ratio 6, and 
a rectangular wing with cambered profile and variable aspect ratio - the 
latter serving as check of the so-called conversion formulas of the air- 
foil theory. 


Nature of the Calculation 

In the aim of explaining the type of calculations employed, the 
simplest case, the straight rectangular monoplane Is to serve as model 
problem. 

The wing is replaced by a continuous surface distribution of lifting 
vortices parallel to the wing center line. On the yawed wing these vor- 
tex lines are diagonal; on the wing with sweepback they are broken at the 

1 "Der Eindecker als tragende Wirbelf lache . " Z.f.a.M.M., vol. 5, 
no. 1, Feb. 1925, PP- 36- 47. 
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apex. The deviations of the wing- profile from the x - y plane are to 
he so small that - as in Birribaum’s report - the lifting vortices in the 
profile may he replaced by those in the x - y plane (the x-axie in 
wing center line toward the right, the y-axis in flight direction 
rearward, and the z-axis at right angles to both downward being counted 
positive) . This assumption is identical with the assumption of Prondtl's 
airfoil theory, according to which interference velocities are small 
relative to the flow velocity. Quantity j(x, y) is taken an the cir- 
culation per unit chord at the point (x, y, 0); that is 

t/'s 

r(x) « // (x, yj dy (l) 

-t/a 

is the circulation distribution over span b, if t is the chord of the 
rectangular wing. To simplify the calculations, the reduced coordi- 

nates \ = —7 — and tj = -k— are used, so that 
b/3 t / a 


+ 1 

r(| ) - J y( t , n) d') 

- 1 

where 7 ( i , tj) - t/27 (x, y) has the dimension of a circulation. For 
7 ( t ; 1 ) the ellipitical distribution over the span 


r ( ti) - y 0 (ti) Ji ~ f 2 


xvas chosen, and for the chordwise distribution the same functions used 
by Bimbaum in his treatment of the two-dimensional problem, 


(a) 7 0 (t) ) = a c 


Fi 


1 + T] 


(b) 7 0 (h) - *> 0 A - 


(c) 7 0 ( n) 


G o r ! J 


According to Prandtl's theory, free vortex lines of intensity 


G 


( i ) 


Sr 

d| 


pass from the lifting vortex surface rearward in flight 
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direction. The downward -velocity w( £ , t]) is computed from y( | ; rj ) 
and e ( | ) by the Biot-Savart law, and from it the profile of the 
wing by the formula 


z 






(2) 


where V is flow velocity. Lastly the drag is obtained by the formula 



( 3 ) 


For the chordwise distribution by this formula the drag then gives the 
suction force S, which occurs at the leading edge of the wing; the 
drag itself is obtained by putting 


7 0 (ti) ° 


1 - ti 

1 + TJ + tt 


and then passing from a to the boundary 0. 

To obtain v( \ , rj) the downward velocity 


w = 


+b /s «> N . 

1 p p oy (x - x' jdx'dy 

*** J J Ox' y(x - x' ) 2 + (y - y r )' a 3 
-b /a o 


kn 


+b /a 

/ 7(x') 


y' dx' 


J (x - x' ) 2 + y' 2 


( 4 ) 


-b /a 


is computed for one lifting filament placed in the x-axis (fig-.,- 2) 


i 1 1 u 
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at point A (with th© coordinates x* and y 1 ) behind the filament. 

The double integral stems from the shedding vortex system (free vort:! ces) , 
the single integral fx'om the lifting vortex itself (bound vortices) . The 
integration with respect to y followed by the introduction of £ and 
tj . gives 

+ 1 + 1 - 

_ 1 [~ 4 c3 1 d £ 1 v t • C>7 d £ ' 

2itb J s i* t - | * + 2«b 2 f si* ( e _ %')j{ 1 _ 

-1 - 1 

i- 7)1 1 P \ 7( £ *' )d V 

2l tb 2 j y ( + n ! 2 ts/b 2 3 

- X 

and after entering of the elliptic distribution y - 1 Q J^- ' £ 2 : 

7 o ^’ t7 o r 1 i • & i • 

w = — + / — — mrrrrr: — - 

2b 2*b 2 j ( | _ £•),/!_ e ,2 </( e - t ') s + ii ’ 2 t 2 /b 2 

- 1 

n't 7 0 ,i‘ .y~Y - V‘ d e 1 

2 ” !,s •/ /7T- e ') 2 + n' s t 2 A a * 


For -rj * = V there is obtained w = 7 0 /' 2 b ('formula (22a) of Prandtl 1 s 

theory). The two additional cylinders result in elliptic integrals, the 
calculation of which is avoided by a special type of series development. 
The task is restricted to the central part of the wing, where 

yr- £ 2 » tj 1 t/b, which obviously is so much better complied with, a 3 
the aspect ratio b/t is greater. The choice for the point A is an 
area ( £ - 5, £ + 5) , such that 


J 1 - W » 6» tj 1 t/b 


( 5 ) 


Within the range ( £ - 5, £ + 5) quantities y and — — are developed 

d£ 
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from | in a Taylor series; outside the range where 

(v t/h) 2 «( e-r) 2 

the roots in the denominator are developed binomially. Assuming uniform 
convergence of the series the integration is made, anticipating that in 
the correlation of the integrals with respect to the inner and outer re- 
gion the quantity 5 cancels. The series are broken off behind the 
squared terms, so that the terms progressing with 

^ f n = 7, 9, ... in the first integral \ 

J 1~- F 2 K n = 5, 1, ... in the second integral'' 


in the integrals with respect to the inner range cancel. As a result, 

all terms with > — : — h from n = 7, and n = 5 are omitted in the 

respective integrals with respect to the outer range. In the binomial 
series all terms with (t)' t/b) n , from n = 3 on, drop out. Therefore 
the true values of the respective integrals with respect to the inner 
range up to the squared terms in q 1 t/b are developed again. The re- 
sult Is: 
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After proving that 6 actually drops out, the squared terms in 
T) ' t/b are disregarded also, which leaves only the linear terms. 


r B 0 l 

w = 7 0 j^A 0 +~ + C 0 T) ' + E q t] ’ In n ' 2 j (6) 


with 


iihmii ii ii ibiiii n^n i 


iii mi iniii i ii i mini i uni i mi nm i nm 
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Substitution 

respect to rj 
inserted for 

formulas for 
(14) from the 


A o = 2 b *■'' 


B ° 1 


o- — | 2 + In J 

27tb -/ 1 - e 2 \ 2 4b(l - 1 2)/ 


-= 


V (7) 


E o= - 


hjtb 2 


-A - E 


.y 


of ti - rj * for n' in (6) followed by integration with. 

' ' from -1 to +1, with one of the three basic functions 
y 0 gives w(| , t]) for the entire superficial wing . The 

induced velocity and profile (the new coefficients derived by 
old ones) are as follows; 


(a) w(| , n) * a [A ° + B°ti + C a 0 T} 2 ] 


z =•■ - 


A a°^ + - ® a V + - C aV 


(b) w(e , n) = b 0 [a^ 0 + B b ° n + i^V] 


z = 


Vi + \ Vi 3 + l Vi 1 


v 

> 


( 3 ) 


(c) w( i , Tj) = c Q [A c ° + C c V+ E c °n 4 ] 


z 



C„°T1 3 



E C°B 5 


J 
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The calculations for the yawed and the swept-back wing are made in 
the same way. The calculation for the wing in sweepback setting must be 
made, aside from (5), with the limitation (2£ tan P + q ' t/b) | - | ') > 

where £ > 0, £' * < 0, and p angle of sweepback* as a result, the .re 

are no data for the middle of the wing. 


Results of the Calculations 
In place of (6) there is in each case an equation: 


w = Jr 


A + — + Cn * * B In r ] 1 2 + Etj * ' in t> 

n 


( 9 ) 


The coefficients are different in each case; those for the straight rec- 
tangular wing carry the subscript o (cf . (7) ), those of the yawed wing 
a cross bar, those for the swept-back wing an asterisk, and those for the 
slipping rectangular wing no subscript. Since the angle P of yaw and 
sweepback is small, everything is developed in powers of (3: 


A = 


2rtb 


4b(l - £ p ) 




rt - 2p ( sin 1 £ In 


B ._n- 


c = 


2rtb a 


J 1 - V 


1- | « t in mil 
.2 t 


> do) 


4b (l - | 2 ) \ 
+ P 2 ( 1 + t s - In ) 


D = - 


V i _ -|2 


p, E 


4rtb 2 | 2 


(1 - P 2 ) 


j 
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For (3 = 0 equations (10), (ll), and (12) change to (7). 

For the subsequent integration with respect to rj * (see above) the 
same formulas are obtained in each case : namely, 

(a) w(| , ti) = a Q [A a + B a ri + C^ 2 ] 

00 v(| , ti) » b D [A b + B bn + C^q 2 + I^q 3 ] (13) 

(c) w( i , q) = C o [A c + B c q + C c q 2 + I^q 3 + E^ 4 ] 

The coefficients with the subscripts a, b ; and c derive from the 
others by the formulas: 
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} (a) 


f V = * 


B a - * 


A + B + - C + D In ^ + E^- + In 

2 4 V2 2 J 


Q + 2D + E( 2 + In ^ , C a = rt E 


Av = % - A + DifiB - - - ^ 
D 2 \ 2 2 / 


1 /J 

(b) ( B i = lt b + - c + Ef- + aift 

i L 2 'S 


1\ 
2 ) 


C’o = « B, 


1A = n - E 

D 3 


” (14) 


A c = 11 


11 
L i 8 


1 \ 


- B -- C + e(- 3m 2 - — ) 
44 l6 y 


r i n 

(c) J B c = n [- D ] C” = « I B - - E 


L‘ 


B r = it - D 


E„ = it t E 

c 6 


From (13) the formulas for the profile z follow immediately. The drag 
is computed by ( 3 ); if A is the lift and q. = - pV 2 the dynamic pres- 

A 2 ^ 

sure, the formula in all cases reads ¥ = — , although the formulas 

rcqh 2 

for w are no longer valid at the edge. The result was to he expected 
by Monk’s stagger law, according to which a variation in the stagger of 
the lift elements is without effect on the drag provided the lift dis- 
tribution of the individual parts is not changed thereby. . 


The same calculations were carried out for the slipping rectangular 

wing with the circulation distributions £ n </ 1 - | 2 (n - 1, 2., 3). 

over the span. The coefficients (in correspondence with n) carry an 
upper subscript. 

For simplification the following abbreviations are used for the 
distributions and their derivatives: 
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* x - t y t 


- f « _ 1 - 2 I 2 


f - _ g e 5 - 3 


fa * -A - I 2 f 2 ’ =1LL1L 

yi - 1 2 


„ 2-9£ 2 + 6j; 4 

Ip == — ■ ■ — 

A- e 2 ' = 


f . = ! 3 / 1 . ?2 f t _ 3 t s - ^ £ 4 

5 yi - 1 2 


f ,! 


1 _ £ 2 3 


Then: 


l (i) 1 j\ 

~ 2to ^ 


- pf hi sin* 1 1 + 2fi’ In - ) 


J ! 


B (l) - 1 f 


C (l) =— f - 2 sin -1 
2to 2 


. fl .. ( i + 1K jhi - I !> -) + _Lihil 


+ P*2 


p 2 ^ - 2 sin" 1 % + fx” C '2 - 1 


tod . i 2 ) N e - 2 e 3 'll 

— t — J *7T=WH 


D ( 1) - -L 3 fl’, 

2rtb 


E (l) * fi”(l + P 2 ) 
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= _L 

2 lib 


«( 3 £2 - p((6-| 2 - 1) sin- 1 ! 


+ 2f 2 ' 15a 


.. 4bd-e 2 ) 


2f 


o 


b< 2 > = 1 f 

nt 


2} 


,(2) _ t 


2rfb 2 l 


6 1 sin' 1 ! 


41(1 - t s ) ^ 3-4|* 
— ) + — + 


. f 2 »( _ + m - s --4 } + — — + p « 6 | 

v 2 t / yi - % s 3 


+ P 2 ; - 6 | sin -1 | + f 2 " (^2 - 


In . 


4b (1 - j 2 ) X 
t y 1 


. 3 - 4 i 2 "1 ] 

+ i > 

7i - V 


2 3 J J 


D^ s) = JL 3f 2 * 
2*b 2 


E<a) =W^ + ^ 
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A (s) = J_ i «(4| 3 - $) - pf(8 I 3 - 2|)ei3 

2rtb [_ V 


4b(l - t 2 ) 

+ 2f ' In i — - - 2f 3 * 


» ±- f 

jtt 15 


’( 3 ) — | - (12 | 2 - l)sin _1 * 

L 


t - t" ( i * in 


41(1 - | 2 ) 


2hb 2 ! 3 " 3 V 2 


t J 


yi - 


r 4l ( l - £ 2 ) \ 

+P 2 ’ - (12 | 2 - l) sin -1 | + f 3 " ^2 - In j 


> (17) 


10 I - 17 i 3 + 6 t l 


D 


( 3 ) _ X 


2jrb 


V 1 - I 


Pf 3 S 


p 2 3 J . 


3 ^ =— - — f 3 "(l + P 2 ) 

4rtb 2 


The calculation gives the following values for lift A, moment 
about the x-axis M X; moment about the y-axis My, drag W, and 


suction force S: 
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7 (t) 

• /l - T] 

J i + Tj 




(b) 



n -/ 1 - n 2 



*ff 2 A pb va d 

it 2 /8 pb 7 b Q 

0 


1 yi -12 

0 

0 

0 

A 

e^yrrp 





« a /l6 pb Va 2 

£ 2 /32 pb Vb 2 

0 


! 3 /i - i 2 

0 

0 

0 


7i - 1 2 


0 

it 2 /64 pbt Vc D 


1 7i - t 2 

0 

0 

0 


lVi - 1 2 





-« 2 /64 pbt Va 2 

0 

it 2 /256 pbt Vc 2 


£ 3 /1 - 1 2 

0 

0 

0 

i 


Vi - e 3 

0 

0 

r ~ ' 

0 

My 

£ 7i - 1 2 

rt 2 /32 pb 2 Va x 

it 2 /64 pb 2 Vb x 

0 


1 2 yi - e 3 

0 

0 

0 


e 3 yi - 1 2 

« 2 /64 pb 2 Va 3 

« 2 /l?8 fibp. Vb 3 

0 


1 - i s 

« 3 /8 p a Q 2 

it 3 // 32 p b Q 2 

0 

W 

i V 1 - £ 2 

it 3 /l6 p a x 2 

n 3 /6Pp b x 2 

0 

iVi- I s 




jc 3 /32 p a 2 2 

it 3 /l 28 p b 2 2 

0 


! 3 7i - I 2 

3 n 3 /l 28 p a 3 2 

3 it 3 / 5 12 p b 3 2 

0 


/TT-pr 

2 Jt/3 b/t p a 0 2 

0 

0 

s 

1 <1 - 1 2 

2 n/15 b/t p a x 2 

0 

0 


1 = 71 - S 2 

2 rt/35 b/t p a 2 2 

0 

0 


| 3 /l- £ 2 

2 n /63 b/t p a 3 2 

0 

0 

i 


The values a 0 , b 0 > o C} a x , . . . b 3 , c 3 are the circulations; 

the letter denotes the chordwise distribution, the subscript the span- 
iri.se distribution. 




















16 


MCA 014 Wo. 1111 


Discrepancies "between the Theory of the Lifting Surface and the 

Theory of the Lifting Line for the Straight Monoplane 

For the drag of the straight rectangular wing "both theories give , 
as already mentioned, the same value. The argument here concerns the 
deviations in the flow or in the profile . For the theory of the lifting 
line also affords a profile of the wing, when the results of the two- 
dimensional problem of the lifting surface are used. The velocity is 
r o 

then — — , increased by the velocity valid for corresponding two- 
2b 

dimensional problem of the lifting surface, the latter elliptical ly dis- 
tributed over the span (reference 1, pp. 20-21, and reference 2, p. 292 ) . 


a_jt r— 

(a) w = + v 1 ~ t 

2b t 


~ a„ 
2 ° 


z = 1 ^2_ 


« + 


V 2 bL t 




r 


(b) w = ~h- + y i - e 2 — n 

4 b t 


z = 


1 b 0 
V" 2 b 


it 2 b 

- T) + 

r\ 1 






(c) w = 0 + J 1 - | 2 — 


l2 - 


1 \ 


1 Ro 


Z = -_ 

V t 


, 1 - 1 N 

2 f , l 3 - - TJ } 

3 2 J 


y (is) 


These formulas agree in first approximation with those of the lift 4* 
ing surface theory; the deviations are of the next higher order in t/b, 
and which are (subscripts F and L referring to values of the lifting 
surface and the lifting line, respectively): 
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X _ p X 

S 3 




(b) z F - z L = 


1 b o t 


8b 2 J 1 - l 2 


r~- 

if- 


3 \ 


a + In 


8b(i - | 2 ) x 


kl9) 


h- 


c n t 


(=, z F - z L - i p±=*(-i+ 

T 8b= yrrp' v *» ^ 


- n 3 -—-'ll® 

3 15 . 




Figure 3 represents the profiles of both theories in the cross sec- 

n 

tions |=0 and | =± The deviations of both curves increase ' 

toward the edge. The extent to which this agrees with reality, or per- 
haps is merely a result of the defective convergence of the series 
toward the edge, is not predictable without further study of this con- 
vergence . 

For the practician, formula (19) affords the following empirical 
formulas for the additional angle of attack Act and the additive curva- 

.1 A 

ture A. £, if c a = is the lift coefficient 

‘ - pV 2 bt 

2 
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Aa = 0.059 - 

h 


A .£ = 0.056 t ° a 


E 


1 t 


Twist of the Yawed Wing and the Swept- back Wing 

From the foregoing it is seen that for elliptic lift distribution 
the wing situated at right angles to the flight direction must receive 
an elliptical spanw3.se distribution in first approximation. If the el- 
liptic lift distribution is then to be retained for yawed or -sweepback 
setting, a new corresponding twist of the wings must be effected. The 
extent of this twist is not immediately obtainable from the given for- 
mulas. From the profiles that were computed for the aspect ratio ^ = 6 

and a number of angles it is seen that the twist increases proportionally 
to the angle 3, for the yawed wing proportional to | and for the 

wing with sweepback approximately proportional to ( - - | 1 1 !. It also 

\ 4 / 

is apparent that the twist must be proportional to the angle of attack 
in first approximation. For on considering a flat plate with zero angle 
of attack it obviously need not be twisted for yawed or sweepback set- 
ting in order to maintain the same flow conditions and hence the ellip- 
tic lift distribution, which in this case is, to be sure, zero. The 
relationship between twist and angle of attack is inconvenient in prac- 
tice, because the twist in a completed airplane cannot be varied with 
the angle of attack. So in a given case the twist must be effected for 
a medium angle of attack. The numerical data for aspect ratio 6, on 
the other hand, afford the following empirical formulas, which should be 
quite useful in practice; that is, for the ordinary profiles for which 


and 


b o > 


0 , 


c o < 


0, and 


and 


are small with respect to 


l o* 


The angle of attack for the wing in yawed setting must be increased 
by about 1.3 iH percent, for the' wing in sweepback setting by 

1.6p fk - I 1 1 ^‘percent of the geometrical angle of - attack (computed from 


the zero lift direction) in the respective cross section without yawed 
or sweepback setting, if the positive course of the yawed wing ( £ > 0)* 
is pushed forward and 3 is the angle, in degrees. The elliptic spanwise 
lift distribution is then retained. 


The ■ Slipping Septangular Monoplane . 

An airplane moving with one corner leading rather than symmetrical 
to itself is said to slip. Consider this same wing as lying symmetri- 
cally in the coordinate system and the free vortices as traveling to the 
side to infinity (fig. 4). For the first, the calculations show little 
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difference from the yawed wing for lateral boundary parallel to flight 
direction. Computing (as under Nature of the Calculation) the velocity, 
w, the integral 


o ' o 

J = 


V d v 


2 « b 


2 rt b 


'_*/ Ji - i * | n! - tanp) 

1 '• D . / 


appears as an essential contribution. 
J = « 


and 


= rt 


' for if + Ti 1 - tanp 
+ I " b 

| + ■n" 5 tanp 


1 - 


< 1 


for 


JTi + V - tanp 2 . lj 


I + n’ - tan P J> 1 


For the far greater part of the wing the first value applies at 
small Pj but the second value must be allowed for in the shaded areas 


of figure 

wholly > 
Hence , if 

W U ; n) 


h. For in this area 


£ + (t) - V) r tan P 
b 


is partly or 


1, if r| 1 runs through the integration range -1 to +1. 

( \ , t}) lies in the shaded area, the portion due to J to 
is: . . . 
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2 Jtbl 


r 

T]+(| + i)b/t cot P i 

<! / ?o 11 dT »* + f 7 o * 


X 


1 


Tj+(| + i)b/t cot P 


t + (tj - ry') - tan P 
b 




J ( t + (t) - V ) - tan P ) 2 - 1 


dr) ' V for i < 0, 


1 J 

2 rtb' 


!ri+( | -x)b/t cot p I + (t| r f]*) ^ tan P 




1 - 


- 1 


/(! + 


!<V 


(•T) - r)' ) - tan p) 2 - 1 
b 


J ' 7 0 n dV l for t > 0 


T)+(| -l)b/t( c6t P 


( 20 ) 


J 


For all bound vortices the shedding free vortices of which still 
cover the point ( | , r s ) (under certain circumstances to be visualized 
as extended toward the negative t) direction) induce on it the constant 

^o . 

downward velocity — . Or simply expressed, if ( t , Tj) lies in the 

2b 

shaded area of figure 4, the w(| , r|) calculated with J - it must be 

supplemented by 


1 

2b 


n 


| + ( r) - T] ' ) - tan P 

b 


/ * 


■J J{i + (n - T)') - tan p) 2 - 1 
rj+(| +i)b/t cot P b 


dt) ' f or | < 0 


V (21) 


n+( I -i)b/t .-cat P 
1 r. i + (h - t' ) ^ tan P 

/ 7 _ 

2b 


:/' * A 


- — ~ dT] ' for | > 0 


| + (tj - T,' )- tan P) 2 - 1 
b 


.7 
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These integrals are elliptic again. Their calculation is omitted since 
the correction obtained through it cannot be applied, because the simple 
formula for w itself becomes invalid at the wing edges. Presumably 
the corners exhibit very high velocities and hence sharp reversals in 
profiles . 

The yawed wing experiences, for lateral boundary parallel to the 
direction of flight, if not twisted conformably to the yawed setting, a 
moment about the y-axis. In figure 5 the half with greater lift is in- 
dicated by +, that with smaller lift by -. The same applies to the 
slipping rectangular wing, as long as the two shaded comers are disre- 
garded. The correction term for the downward velocity, which these 
corners furnish, then becomes negative for % > 0 and positive for 
| < 0; hence, to obtain elliptic lift distribution, the left corner 

must be set steeper, the right corner not so steep. . (perhaps even steep 
downward) . If this is not done, the left corner gives less lift, but 
the right one more lift (fig. 5). Therefore it is impossible to say in 
what sense the moment about the y-axis of a symmetrical rectangular 
wing in yaw occurs, because without calculation it is not known which of 
the two effects we dominates . 


Second Problem for the Flat Monoplane of Aspect Ratio $ 

For the practice the second fundamental problem, namely, the deter- 
mination of a lift distribution for a given profile, will always be the 
more important one . 

This problem is treated by assuming a circulation distribution with 
indeterminate coefficients, calculating the flow, profile, and drag for 
every single distribution as for the solution of the first fundamental 
problem, and then determine the indeterminate coefficients from the con- 
dition that the wing shall assume a prescribed profile. The wing copies 
the profile exactly at as many points as there are indeterminate coeffi- 
cients in the lift distribution. 

For the slipping rectangular flat wing the following circulation 
distribution with 12 coefficients is assumed: 
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= J 1 - % 2 f j a D /,- 1 '- ~- + b Q ,/ 1 - T 1 S + CqT! y 1 - 


. </ 1 + T) 


+ i 


a i 


^ + bi */ 1 - T| 2 + CiT) «/l - 


^ 1 + T] 


iT I 


+ 


&2 


+ b 2 yi - TJ 2 + C 2ri 7l - r f 2 


1 + 1} 


+ d 


l w 


/A-AjJ + t 3 y 1 - q 2 + C 3 T] : 1 T T ] 2 ] [ 

1 + q J J 


The results of the necessary calculations, which correspond to 
those of the first fundamental problem, have already been described un- 
der Results of the Calculations. Velocity, profile, lift, and moment^ 
are linearly composed of the individual parts; the drag is quadratic in 
the circulation and is therefore tedious to compute in certain cases. 
That the wing shall be flat simply implies that the velocity w( \ , q) 
over the entire wing shall be constant. Thus at 12 points of the wing 
the equation: 



W U , Tj) = a 0 [A a + B a T] + C a T ] 2 ] + b Q [A b + B b T] + C ^ 2 + ^ 3 ] + c Q [ A c + B Q r\ + C c T ) 2 + D ^ 3 + E cn 4 ] 


+ ai[A a (l) + + C a * 1 V ] + M A b (l) + ] + Cl [A c ( ^ + ] 

+ a 2 [A a ( 2 ) + ] + b 2 [A b ( 2 ) + ] + c 2 [A c ( 2 ) + ] 

+ St A a (3) + ] + l> 3 [A b (3) + ] + c 3 tA 0 (3) + ] 

=' constant (22) 

must be fulfilled. 


ro 

u> 
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The constant is proportional to the angle of attack a,. The coef- 
ficients of equation (22) follow by means of (14) from (12), (15), (l 6), 
and (IT). The intersection points of the straight line f = ± - , ± - 

, _ 4 4 

with the straight lines rj - 0, ± - Jz are chosen as the 12 points at. 

2 

which equation (22) is to he fulfilled. These values are chosen, be- 
cause any contour is better approximated by one such third order with 
the same tangent angles in the points rj = 0, ± i J ~ 3 than otherwise by 

a curve of the third order (reference 2, p. 294). The solution of the 

- no s 


12 equations for different values of P 
lowing values: 


(*■ 


2 6 V tan a, J 


j affords the f ol- 


p 





mi 



b 2 X 

c 2 X 

£Lrj W 

b 3 X v 

; C 3' 4 

0 


-1.46 

-0.23 

0 

• 

0 

0 

29.5 

-20.5 

7.35 

0 

• 0 

'0 

.1 

66.8 

-1.46 

-.17 

-1.58 

.09 

1.07 

29.4 

-20.1 

6.68 

-1.77 

: 3772 

p.l>95 

.3 

66.7 

-1.14 

-.25 

-4.58 

-1.01 

1.43 

30.5 

-20.1 

7.30 

-5.40 

11.4 

-2.90 

.5 

66.5 

-0.48 

:.29 

-7.54 

-1.76 

1.68 

31.5 

-20.0 

7.37 

-9.02 

19.0 

-3.48 


The spanwise circulation distribution is: 


7 = 66.07 * J\ - t f 1 
= 66.07 tt „/ l - i 3 [ 1 

= 66.13 * y i - r r i 

= 66.21 « V 1 - I 2 [ 1 


0.023 i 
0.077 t 
0.127 £ 


+ 07292 t ] for p 
+ 0.293 I 2 + 0.0014 £3] • for p 
+ 0.309 + 0.0045 I 3 ] for P 
+ 0.325 i 2 + O.OO76 | 3 J for P 


0 

0.1 

0.3 

0.5 


(See fig. 6.) 


The coefficient of i s for P = 0, according to Bets (reference 
3), for a rectangular wing of identical profile and angle of attack is 
in round figures O.32O (found by graphical interpolation for L = 3.82 
corresponding to b:t = 6), which is in sufficient agreement with the 
value obtained here. The discrepancy is due to the fact that the series 
was broken off behind the term with | 2 f whereas Betz completed the 
whole series. 


The moment about the y-.axis introduced by the yawed setting is: 
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My = - — p b 3 T 2 tan a X 0 

for (3 = 0j , 

J 16 


x 1.50 x icr 3 

for 0 = 0.1 

X 4.93 X 10- 3 

for {3 = 0.3 ; 

x 8.17 X 10- 3 

for (3 = 0.5 


My increases proportional to 3; the greater lift is experienced by the 

trailing half of the wing. However, it should be noted that a test may 
well give the opposite result, because the two shaded corners (fig. 4) 
have been disregarded; they presumably have, a considerable share in the 
formation of the moment. 

The calculation of the total drag (with suction force) is consider- 
ably reduced for the wing in question, since w « V tan a is constantly 
distributed over the span and chord of the wing. It simply gives 


¥ + S = A - 
V 


= A tan a . p 

2 


b 2 V 2 tan 2 a x 
x 
x 
x 


70.88 x io~ 3 
70.91 x 10“ 3 
71.24 x 10“ 3 
71.59 x 10“ 3 


for 

P = 0 

for 

0 = 0.1 

for 

0 = 0.3 

for 

P = 0.5 


The separation in drag and suction force is effected by computing 
the suction force alone; only the chord distribution (a) in equation 
(3) and the term of the speed arising from B are considered. 


p b 2 V s tan 2 a X 1 6.6 x 10 -3 

for 

P 

= 0 

x 16.7 X 10“ 3 

for 

P 

= 0.1 

X 16.7 x 10~ 3 

for 

P 

<= 0.3 

x 16.8 x 10- 3 

for 

P 

= 0.5 
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For $=0 the theory of the lifting line (reference 1, formula 
(23)) gives: 


W 




- ( a.? n + b 2 ( =— p b 2 V 2 tan 2 a x 16.CO x 10 3 

8 V 2 2 2 S J 2 


in satisfactory agreement with the value computed here. 


Conversion Formulas 

The theory of the lifting line gives the formulas (reference 2, 
formulas (26a) and (27a)): 


o>i - - 


ca 


V b-, 



C-Wl - CV2 


It. (h 

« \ hi b 2 /’ 


(23) 


■with Co = 


q b t 


and 


w 


c w = 


q b t 



dynamic pressure. 


These formulas make it possible to convert the data from one case (l) to 
another (2) with different aspect ratio for fixed cc a ; that is, to ap- 
ply the test data obtained for one aspect ratio to other aspect ratios. 
These formulas were checked for a specific profile (Gottingen Ho. 3&9, 
fig. 7) up to b/t = 7 to 1, and found well confirmed up to b/t = 3; 
only those for b/t = 2 and 1 manifested considerable discrepancies 
(reference 4). It is suspected that these discrepancies are confirmed 
by the theory of the lifting surface. 

The profile of the wing (that is, its backbone) is replaced by a 
curve which aerodynamically is as equivalent as possible to the profile. 
As such, a center line as symmetrical as possible is chosen. The calcu- 
lation, reproduced below, did not give the desired results; the diagram 
c a (c w ) exhibited slight discrepancies with respect to the lifting line 



NACA TM No. 1111 


ST 


theory; the diagram c a (a) shoved, on the whole , substantially lover 
lift coefficients than the test. This result indicated that the hack- 
hone was not quite exactly chosen; it should perhaps he taken nearer to 
the suction side, the predominant effect of which is known. 

In order to establish the ratio of the distance of the backbone 
from the upper and lower edge of the profile for the backbone to be 
aerodynamically equivalent to the profile, a number of circular crescent 
Profiles that give equal lift for equal chord were analyzed by the 
Karman and Trefftz method (reference 5) . The looked-for ratio is de- 
pendent upon the peripheral angles of the two circular arcs of the 
crescent. The present profile was related to a' circular crescent, and 
the ratio determined for these profiles was about 2/3. The backbone 
obtained with this value is represented in figure 7 . 

Now the lift distribution with the six coefficients a Q , b Q , c D , 
a 2 > b 2 , c 2 for varying .b/t is computed, as under 6, for a wing with 
this constant profile. To this end the angles of the tangent to the 
backbone are measured for 

■n 



with respect to the chord drawn below the backbone. 


n 



13° >5’ 


q = 0 : eg * 2° 7’ 


q = -JL.: a 3 « 6° 23’ 30" 

2 • 

At the points % = ± - ± 2; q = 0, ± r! — the equation 

4 4 , 2 

w ( I » Tl) = a o [Aa + B a q + C a q 2 ] +. b Q [A’ 0 + Bfc q + q 2 + Db q 3 ] 
+ c q [A c + B c q + C c q2 + D c q 3 + q 4 ] + a 2 [A a (2) + B a (s } q 


( 3 ) 


+ C a V2; T) 2 3 + b 2 [A b ^ 3 ^ + ]+ .< s £ a c (3) + ] = V(cc + a,,) (2^) 


(3) 



must tie fulfilled, with a - 0 indicating one of the three values a 2 , a 3 depending upon 

the value of T) . The solution of the equations gives: 


10 3 

10 3 

■U - 

0 2 t V 

0 2 t V 

- 21 . 8 + 400.8 a 

100.0- 8.76 a 

- 24 . 0 + 360.0 a 

99.6- 22.92 a 

- 28 . 8 + 278.4 a 

96.8- 66.0 a 

-31.8+211.8 a 

.. . i 

90.6-120.8 a 


10 3 

10 3 
V, 

10 3 

dp 

2 t 7 

Dp 

t v 2 

Cp 

2 t V. 

- 32 . 2 + 176.4 a 

85 . 0 - 123.0 a 

-65. !h- 37.7 a. 

- 30 . 6 + 138.4 a 

76 . 4 - 151.2 a 

-67.6+ 74.8 a 

-25.6+ 92.4 a 

56.0-137.6 a 

-68.2+168.0 a 

- 19 . 1 + 34.4 a 

33 . 5 - 23.3 a 

-62.2+210.3 a 
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From this c a and. c w are computed as functions of a and com- 
pared with the test data, while it is home in mind that a is referred 
to the chord below the profile j the difference is 1°9'20". Regarding 
the values obtained here by means of the .lifting surface theory like ex- 
perimental data; that is, converting them by (23) to the aspect ratio 
os, produces no substantial discrepancies in the d a (c v ) diagram from 

the lifting line theory, or, in other words, the quoted conversion for- 
mula for c w is confirmed. The c a (a) diagram exhibits good agreement 
with the test, at least within the range of small angles of attack (fig. 
8). The differences between test on conversion formula for a at 
b/t = 2 and 1 are largely confirmed.. Obviously another effect disre- 
garded so far as. the theory is involved that causes the curvature of the 
curves c a (a) and appears to be of significance especially for small 

aspect ratios. Its explanation with the available means cf the airfo.il 
theory is probably impossible, as the latter is essentially a linear 
theory. 


Translation by J-; Vanier, 
National Advisory Committee 
for Aeronautics. 
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Figs. 6,7,8 



Figure 6.- Lift distribution of the flat rectangular 
wing of aspect ratio 6 at 0° and 28.6° , 
yawed setting. ' 


